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S : ABSTRACT 
(N 

During the last few years the Tevatron has dramatically improved the bounds on rare i?-meson 

00 : n + - 

. decays into two leptons. In the case of — > , the current bound is only ten times greater 

■ than the Standard Model expectation. Sensitivity to this decay is one of the benchmark goals 

^ ' for LHCb performance and physics. The Higgs penguin dominates this rate in the region of large 

H : 

^ tan P of the MSSM. This is not necessarily the case in the region of low tan /?, since box and Z- 

penguin diagrams may contribute at a comparable rate. In this article, we compute the complete 
one- loop MSSM contribution to B'^ ^ ^ i~^£'~ for = e,/i. We study the predictions for general 
values of tan /? with arbitrary flavour mixing parameters. We discuss the possibility of both 
enhancing and suppressing the branching ratios relative to their Standard Model expectations. In 
particular, we find that there are "cancellation regions" in parameter space where the branching 
ratio is suppressed well below the Standard Model expectation, making it effectively invisible to 
the LHC. 



1 Introduction 



One of the most promising signals for new physics at the LHC is the rare decay — > jji^ jjL~ . 
This decay is suppressed as a loop-level flavour-changing neutral current and by a lepton 
mass insertion required for the final state muon helicities. The LHC will be the first 
experiment to be able to probe this decay channel all the way down to its Standard Model 
(SM) branching ratio. The decay is especially 'clean' because its final state is easily tagged 
and its only hadronic uncertainties come from the hadronic decay constant fs^- Further, 
enhancements to this branching ratio by new physics can be resolved with only a few 
inverse femtobarns of data, making this an exciting channel for beyond the standard model 
searches in the first few years of LHC operation. 

The current experimental status and the Standard Model predictions for the branching 
ratios B{B'^^ to leading order in QCD are displayed in Table [H This is the 

updated version of the Table 1 presented in the review of Ref. [1]. Further reviews can be 
found in Ref. [2]. 



Channel 


Expt. 


Bound (90% CL) 


SM Prediction 




CDF II [3] 


< 4.7 X 10^* 


(4.8 ± 1.3) X 10-^ 




CDF II [3] 


< 1.5 X 10-8 


(1.4 ±0.4) X 10-10 




CDF [4] 


< 6.1 X 10-6 


^ 




BABAR [5] 


< 9.2 X 10-8 


^ 



Table 1: Current experimental bounds and SM expectations for leptonic i^^-meson decays. 

The Standard Model predictions for the dimuon decay of B^ and 5° mesons were first 
calculated by Buchalla and Buras in Ref. [6] and Higgs penguin contributions in Ref. [7]. 
Their analysis can be generalised to include lepton flavour-violating decays with the final 
state n~^e~ which are not measurable within the SM extended with see-saw neutrino masses. 
The error in the SM predictions for the 5^^^ branching ratios originates primarily from the 
uncertainties in the decay constants [8], 

fs, = 230 ± 30 MeV , /b, = 200 ± 30 MeV , (1.1) 

and in the top-strange and top-down elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [9], 

I l^t, I = 0.0406 ±0.0027, = 0.0074 ± 0.0008 . (1.2) 
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The dimuon decay of is of particular interest to experimentalists because it is a 
benchmark process for LHCb physics and performance. The LHCb will be able to directly 
probe the SM predictions for this rare decay mode at 3a {5a) significance with 2 fb~^ (6 
fb~^) of data, or after about one year (three years) of design luminosity [10]. In addition, the 
general purpose detectors ATLAS and CMS will also be able to reconstruct the — > //''"//~ 
signal with significance of 3a after Ri 30 fb~^ [11]. It is not clear whether LHC can reach 
the SM expectation for B^ — > jjf^jjr . 

At the dawn of the LHC era, it is important to understand the possible contributions of 
new physics to a discovery in the B^^ — > channels. These are particularly promising 
decay channels for the Minimal Supersymmetric Standard Model (MSSM). Under the as- 
sumption of large values of tan/3 and Minimal Flavour Violation (MFV), where the CKM 
matrix is the only source of CP and flavour violation, the branching ratio for B^ — > /z"^//"" 
is dominated by the Higgs penguin mode and is approximately given by 

v-).wo-'(^|i^y(55i|fX)', ,,3, 

where Ma is the CP-odd Higgs mass. Thus in the large tan/3 regime this branching 
ratio can be significantly enhanced over the Standard Model expectation. This has been 
discussed extensively in Refs. [1,12-15]. The large tan /3 regime is preferred, for example, by 
supersymmetric 5'O(10) grand unified models. Further, the currently observed excess in the 
anomalous magnetic moment of the muon (gr— 2)^ [16] implies an additional enhancement in 
B{B^g — > jjf^jjr) in certain supergravity scenarios [17]. A field theoretic study of this decay 
in the large tan (5 limit focusing on the resummation of tan (5 was conducted in Ref. [18-21]. 

Thus far, however, the published analyses have focused primarily on the large tan (3 
region with MFV and have neglected possible flavour mixing in the squark sector. With 
the upcoming experimental probes of -6(5°^ — > t^l'') down to the SM expectation, it 
is important to undertake a full, general calculation of this branching ratio without a 
priori assumptions on the pattern of squark and slepton flavour mixing or electroweak 
symmetry breaking. In particular, in the region of low tan /3, the effects of box and Z- 
penguin diagrams could be of the same order as the tan /3-enhanced Higgs penguins. The 
interference of these terms could conceivably lead to a cancellation that would suppress 
the branching ratio helow the SM prediction. This region of parameter space has not yet 
been thoroughly investigated. This paper fills the gap in the hterature on the low tan /3 
properties of these decay modes. 

If the branching ratio is significantly enhanced by new physics, it may even be visible 
at the Tevatron. Alternately, if it is significantly suppressed by new physics, it may be 
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mvisible even at the LHC Either way, the status of this decay could become an important 

factor for planned LHCb upgrades. For example, it could play a critical role in determining 

whether an LHCb upgrade should focus on a more precise measurement of 

or instead reach for the branching ratio of /^^/^~ which is an order of magnitude 

smaller. 

In this article we calculate MSSM predictions for dileptonic i?^, B^ decays with arbitrary 
flavour mixing. In our numerical analysis, we ignore r-lepton final states since decays like 
B^ T^T^ or 5° r+yU~ since they cannot be observed accurately at the Tevatron or 
LHC. Although our calculation is sufficiently general to include lepton flavour violating 
5-decays like B^ /W^e^, we do not consider them in our numerical analysis due to their 
small branching ratio ( < 10^^^) at low tan/^. We therefore concentrate on the decays 

— > fi^fi^ and B^ The general calculation is presented in the appendix and 

the code used in our numerical analysis is available to the public^. 



2 Effective Operators and Branching Ratios 

There are ten effective operators governing the dynamics of the quarks-to-leptons transition 
qiqj _^ £+K £-L^ ^j^]^ ^ ^^q2 ^ g^q3 ^ ^^^^ £i ^ £2 ^ ^^£3 ^ ^_ ^^le effective 

Hamiltonian reads: 

'H = . . ( CvXY OvXY + CsXY OsxY + Ctx Otx \ , (2.1) 

where flavour and colour indices have been suppressed for brevity. The (V)ector, (S)calar 
and (T)ensor operators are respectively given by 

O'vfy'^ = (^l^Pxq'Wl.PY^''), 
0¥xY = (FPxg')(FPy£^), 

O'tW^ = (^a^''Pxq'){¥cj,,i^). (2.2) 

We follow the PDG conventions for the quark content of the i?°-mesons, 5^ = 6s and 
5^ = hd [23]. Thus in Eq. ((221) we identify q-^ = h and = s oi d for 5°^ f+f" 
respectively. 

The explicit forms of the Wilson coefficients for the MSSM are calculated at the elec- 
troweak scale, Q = Mw- These are given in the appendix. The contributions to these 

^Predictions of MSSM for B^^t or B^fie a.t large tan/3 have been investigated in Ref. [22]. 
^In order to obtain the Fortran code, please send e-mail to jctnusz.rosiek@fuw.edu.pl 
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coefficients can be classified into Z-penguins, Higgs-penguins and box diagrams, shown 
in Fig. [H The photon penguin contribution 5° —>■ vanishes in matrix element cal- 

culations due to the Ward identity. We do not consider the very large tan/? scenario 
(tan/9 > 30), since in this region our calculation has nothing to add to the current litera- 
ture (see previous section for references). Thus no resummation of higher orders in tan/3 
is necessary and all formulae given in the appendix are strictly one-loop. 




We now focus on the decay — > Corresponding formulae for the decays 

can be derived analogously. Down quark vector and scalar currents hadronize to i?s-mesons 
as 

{0\h,PLiR)s\Bsip)) = -(+)'-pJs^, (2.3) 

A Till) + Tflg 

where = pi^ + p^^ is the momentum of the decaying i?s-meson of mass M^^. In 
deriving Eq. (12. 4p we have used the quark equations of motion. One immediate conse- 
quence of the -B° — i> kinematics is that tensor operators vanish in the matrix element 
{0\b(Tfj,u s\Bs{p)) because there is no way to make an antisymmetric tensor with the single 
available momentum p^. This reduces the total number of effective operators contributing 
to Bg — i> in Eq. (12.21) from ten to eight. The matrix element is therefore 

M = Fsle + Fp'hd + Fvp^h/ + FAP^hf^l^i, (2.5) 

where the £s correspond to external lepton spinors, e.g. FsJi = Fsu{p2) v{p4), etc. The 
momenta are assigned as in Fig. [TJ The (S)calar, (P)seudoscalar, (V)ector and (A)xial- 
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vector form factors in Eq. (12.51) are given by 

Fs = -7 — ^'"^ {CsLL + CsLR — CsRR — Csrl) , (2.6) 
4: mb + iris 

Fp = - — ^""^ {—CsLL + CsLR — CsRR + Csrl) , (2.7) 
4 rrif, + rris 

1 

Fv = —-^fBs{CvLL + CvLR — CvRR — CvRL), (2.8) 

Fa = ~^fBs {~CvLL + CvLR — CvRR + Cvrl) ■ (2.9) 

It is now straightforward to square the matrix element in Eq. (12.51) . and determine the 
branching ratio for the decay 

B{B.^i,i,) - ) f-[ ) , (2.10) 

where tb^ is the hfetime of Bg meson, and 

\M\' = 2\Fs\' [Ml - {me, + me^Y] + 2\Fp\^ [M^ - {me, - me,)'] 
+ 2\Fv\' [MlXme, - me J' - {m^ - m\f\ 
+ 2\Fa\' [MlSme, + m.J^ - (m,\. - m^J^] 
+ ARe{FsF*){me, - me,) [M^ + {me, + me,)'] 

+ 4Re{FpFl){me, + me,) [M^ - {me, - me,)'] . (2.11) 

Notice that the contribution from the vector amphtude, Fy, vanishes in the lepton 
flavour conserving case, L = K. In this case the formula in Eq. (12.111) agrees with results 
of Ref. [24]. 

The form factors in Eqs. (I2.6H2.9]) do not receive additional renormalisation due to 
QCD corrections. The conservation of axial-vector current operators {Ovxy) result in 
vanishing anomalous dimension associated with this operator. The scalar operators {Osxy) 
renormalise like a quark mass parameter and thus the ratio Csxy{Q) /[^b{Q) + ^s{Q)] is 
a renormalisation group invariant quantity [25]. 

Wilson coefficients and parameters entering Eqs. (I2.6H2.9]) and (12.111) are all calculated 
at the top quark mass scale, i.e. Q = mt- The quark pole masses, mb or m^, are related 
to their Di?-running one-loop quark masses at the scale Q, mq{Q), by the well-known 
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formulae: 



TUbiQ) 

mt{mt) 



TUf, 



nit 



371 



asirrih) 



bo 



(2.12) 
(2.13) 



with bo = 11 — 2?2j/3 with nj = 5. Since our calculation for the SUSY corrections is 
performed in the DR renormalization scheme [26], our initial conditions for parameters 
must be converted into this scheme. Eqs. (12.121) and (I2.13P contain the appropriate MS — > 
DR conversion factors [27]. Similar conversions for gauge couplings is small and is ignored 
in our numerical results. 

We have included the general decays of Eq. (I2.10p in our numerical code. Due to 
small branching ratios, however, it is unlikely that we will observe lepton flavour violating 
5-meson decays at moderate or small values of tan/3 and so we will not consider these 
processes in the remainder of this paper. 



3 Numerical Analysis of Bs^d 

3.1 Structure of the MSSM contributions 

We now focus on the lepton flavour- conserving processes B^j,—^ fi^fi~. Recall from Section 
[T]that we are interested in cases where the branching ratios for these processes are either 
enhanced or suppressed significantly relative to their Standard Model predictions. An 
enhancement would either lead to an early discovery at the LHC (or even the Tevatron) or 
stronger constraints on the allowed magnitude of squark flavour violation. A suppression, 
on the other hand, could lead to a non-observation of fJ'^fJ'~ at the LHCb due to 

cancellations from new physics in the decay amplitude. 

For the lepton flavour conserving decays £x = = and the squared amplitude in 
Eq. (12. lip takes the form 

\M\''^2Ml{\Fs\^ + |Fp + 2m^F^p) , (3.1) 

where we have also taken the limit m^/MB^ 0. We may distinguish two possible scenarios 
for the relative size of the MSSM contributions to the right-hand side of Eq. (13. ip : 
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1. Higgs penguin domination or large tan (3 > 10 

In this large tan/? regime one can usually expect an enhancement of the branching 
ratios as in Eq. (11. 3p . This case has been thoroughly investigated in the literature, 
although mostly in the limit of minimal flavour violation and vanishing intergenera- 
tional squark mixing. In such case it turns out that I-F5I ~ \Fp\ ^ 2mi\Fyi\ because 
of tan^ j3 enhancements. Although this is the standard situation for large tan (3, it is 
not general since a kind of Glashow -Iliopoulos -Maiani (GIM) cancellation mecha- 
nism may result in F§^^ ^ [18,21], thus making the box and Z-penguin diagrams 
phenomenologically relevant. 

2. Comparable Box, Z-penguin and Higgs penguin contributions or low tan (3 < 10 

In this low tan j3 case the supersymmetric Higgs-mediated form factors Fs^p are sup- 
pressed and become comparable to or even smaller than Fa- Thus the full one-loop 
corrections to the amplitude are needed. These are presented in the appendix. In 
this case either an enhancement or a suppression of the branching ratios is possible 
depending on the particular choice of MSSM parameters. 

Barring accidental cancellations, an enhancement of the branching ratios can come from 
any of the contributions in Fig. [TJ On the other hand, it is a bit trickier to suppress the 
branching ratios below their Standard Model predictions as this requires a cancellation 
between various terms. This is the case we would like to investigate further. 

We would like to find the minima of B{B^^ fi~^fi~), i.e. the minima of Eq. fl3.ip . We 
distinguish between two cases: 

Fp + 2 F4 ^ and Fp > , (3.2) 
or 

\Fs\ ^ \Fp\ ^ IF4I ^ . (3.3) 

In the first case, Eq. (13.21) . the pseudoscalar and axial contributions cancel while the 
scalar contribution is negligible. This can be realized, for example, in models where the 
MSSM is extended with an additional, light, CP-odd Higgs boson. Ref. [28] shows that this 
can occur even in the minimal flavour violating limit of such a model. Such cancellations, 
however, can also take place in the general MSSM when left- and right-handed squarks mix 
in the strange and charm sectors. Furthrmore, it has been pointed out in Ref. [29], that 
interference between the scalar /pseudoscalar new physics and Standard Model operators 
can decrease the i3(-Bs — > /U+/i^) far below its SM prediction. This is explored further 
within MSSM in the numerical analysis of Section [3^ 
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The second case, Eq. (13.31) . happens when Higgs contributions are neghgible compared 
to the axial contribution (i.e. low tan/5 and large Ma) and Fa becomes small due to 
cancellations among the CyxY coefficients in Eq. (12. 8p . Our numerical analysis shows 
that such a cancellation is possible but requires a certain amount of fine tuning once 
constraints on squark mass insertions from other fiavour-changing neutral current (FCNC) 
measurements are imposed. 

3.2 Numerical setup 

To quantitatively study the effects mentioned in the previous section, we perform a scan 
over the MSSM parameter space. The ranges of variation over MSSM parameters are shown 
in Table [21 Because our numerical analysis is based on the general calculation presented 
in the previous section, we are not restricted to particular values of tan/5 or the MFV 
scenario. Flavour violation is parameterised by the "mass insertions", defined as in [30,31], 

5{f,y = (BM=. (3.4) 

As before, /, J denote quark fiavours, X, Y denote superfield chirality, and Q indicates 
either the up or down quark superfield sector. 

To realistically estimate the allowed range for B{B^^ —>■ one must account the 

experimental constraints from measurements of many other rare decays. SUSY mass inser- 
tions, in particular, are strongly constrained by such measurements. The most important 
constraints have been calculated in the framework of the general MSSM using a standard 
set of conventions [19,31-36]. We have used the library of numerical codes developed in 
those studies to bound the MSSM parameter space based on the set of observables listed in 
Table [31 no further bounds (e.g. dark matter, electroweak observables, etc.) are imposed 
other than those listed. 

For all the quantities in Table [3] for which the experimental result and its error are 
known, we require 

For the quantities for which only the upper bound is known, we require 

+ <Q'"^. (3.6) 

The first and second terms on the right-hand side of Eq. (13.51) represent the 3a ex- 
perimental error and the theoretical error respectively. The latter differs from quantity 
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Parameter 


Symbol 


Min 


Max 


Step 


Ratio of Higgs vevs 


tan P 


2 


30 


varied 


CKM phase 


7 





TT 


7r/25 


CP-odd Higgs mass 


Ma 


100 


500 


200 


SUSY Higgs mixing 




-450 


450 


300 


SU (2) gaugino mass 


M2 


100 


500 


200 


Gluino mass 


Ms 


3M2 


3M2 





SUSY scale 


MsvSY 


500 


1000 


500 


Slepton Masses 


Mi 


^susy/3 


^susy/3 





Left top squark mass 




200 


500 


300 


Right bottom squark mass 


Mr 
Or 


200 


500 


300 


Right top squark mass 


Mr 


150 


300 


150 


Mass insertion 


^13 ^23 
"dLL' "dLL 


-1 


1 


1/10 


Mass insertion 


SIS ^23 
"dLW "dLR 


-0.1 


0.1 


1/100 



Table 2: The range of input parameters for the numerical scan. "SUSY scale" refers to the 
common mass parameter for the first two squark generations. The parameter tan j3 takes on 
values within the set: tan/5 = (2, 4, 6, 8, 10, 13, 16, 19, 22, 25, 30). All mass parameters are 
in GeV. The top quark pole mass have been taken 171 GeV. Imaginary part of parameters 
^dLL^ ^dLRy M2 have been set to zero. The trilinear soft SUSY breaking couplings are 

set to At = Af, = Mq^ and Af = Mi throughout. 

to quantity and is usually smaller than the value q = 50% which we assume generically 
in all calculations. Apart from the theoretical errors that come from uncertainties in the 
QCD evolution and hadronic matrix elements, one must also take into account the limited 
density of a numerical scan. In principle, with a very dense scan and sufficient computing 
time, it should be possible to find SUSY parameters that fulfill Eq. (13. 5p within the calcu- 
lation's "true" theoretical errors. This, however, may not be necessary and may even be 
undesirable. Our goal is to find "generic" values for the branching ratio B{B'^^ fi^fi^), 
i.e. values allowed by fairly wide ranges of SUSY parameters without strong fine tuning or 
the need to resort to special points in parameter space where "miraculous" cancellations 
evade experimental bounds. In our scan we thus use wide "theoretical" errors assuming 
that this procedure faithfully represents the ranges of the MSSM parameters. If necessary 
the exact values of parameters fulfilling the bound in Eq. (13.51) with smaller q can be found. 
A more detailed discussion of the problems associated with scanning over multidimensional 
MSSM parameter space can be found in [36]. 
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Quantity 


Current Measurement 


Experimental Error 


mo 


> 46 GeV 






> 94 GeV 




mi 


> 89 GeV 




mi 


> 95.7 GeV 




rrih 


> 92.8 GeV 




|ex| 


2.232 ■ 10-3 


0.007-10-3 




3.483 ■ 10-15 


0.006 ■ 10-15 


|AMb| 


< 0.46 ■ 10-13 






3.337-10-13 GeV 


0.033 ■ 10-13 GeV 


AMb, 


116.96 ■ 10-13 GeV 


0.79 • 10-13 GeV 


Br{B ^ X,7) 


3.34 ■ 10-^ 


0.38 • 10-^ 


Br(iri ^ ttVz/) 


< 1.5 ■ 10-1° 




Br(ir+ ^ TT+z/z/) 


1.5 ■ 10-1° 


1.3 ■ 10-1° 


Electron EDM 


< 0.07-10-26 




Neutron EDM 


< 0.63 ■ 10-25 





Table 3: Constraints used in the scan over MSSM parameters. LEP data are used for the 
Higgs mass bound [37], i.e. mt > 92.8 — 114 GeV depending on the value of sin2(a — p). 

3.3 Predictions for B{Bs A^^A^ ) 

Fig. [2] shows the predictions for B{B^ f^^fJ' ) over a general scan of 20 million points in 
parameter space according to Table [3] and including the bounds described in the previous 
section. The upper bound set by CDF in Table [H depicted as a solid red line, can be 
attained even with very low values of tan (3. We focus on the lower limit of the branching 
ratio and therefore restrict to the region of parameter space where tan/5 < 30. In this way 
we also avoid the technical complications connected with the resummation of higher order 
terms, discussed in Refs. [18-21]. We vary ^J^^^ (upper panel) and ^^^.ij (lower panel) one 
at a time while setting the other to zero, e.g. all ~ ^ ^^"^ only 6"^^^ 7^ in the upper 
panel. 

When S'^\i^ is varied in the range [—1, 1], we find B{B^ f^'^fJ'~)min ~ 10-^. This min- 
imum is almost independent of tan f3 but depends on the magnitude of the mass insertion 
(upper right panel). I^^'^^l can take on values up to ~ 0.9 and still pass all the constraints 
in Table [3l though points beyond 0.3 are less dense. We note here the importance of cor- 
rectly incorporating the LEP Higgs mass bound. If for example we set mh > 114 GeV 
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independently of the value of the ZZH coupling, then is restricted to values smaller 

than ^ 0.3. 
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Figure 2: Upper panel: Predictions for B{Bs versus tan/5 (left) and 5^^^ (right) 

from the scan of MSSM parameters in Table [2] and under the constraints of Table [31 The 
upper solid line shows the current upper bound from the Tevatron and the lower dashed 
line the SM expectation. Lower panel: Similar to the upper panel but with 5^'^^ varied. 



More interesting is the case when 5^"^^ is varied in the range [—0.1, 0.1]. We find a nar- 
row cancellation region around ~ —0.01 and tan/5 < 10 where -23(5° fJ'~^fi~)min ~ 
10~^^ (lower right panel). This is three orders of magnitude lower than the Standard Model 
prediction, making it effectively unobservable at the LHC In order to better understand 
cancellation region we study a representative point with a very low branching ratio, for 
example: 

tan/5 = 4. Ma = 300, /i = -450, Ma = 100, M3 = 300, 
SUSY scale = 400, M^^ = 150, = M^^ = M^^^ = 600 , (3.7) 

where all masses are in GeV. 
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Figure 3: Contributions to B{B^g^ A*^/^") from various parts with the parameters in 
Eq. fl3.7l) . Left: Contributions from the diagrams in Fig. [T] versus Right: Magnitude 

of the form factors appearing in Eqs. fl2.6H2.9l) versus 



The cancellation is easy to understand if one independently considers the contributions 
to the branching ratio from each diagram, as shown on the left in Fig. [31 The 'Box', 'Higgs' 
and 'Z' lines indicate the value of B{B^^ — > f^^t^ ) given by only the listed contribution with 
all others set to zero. The total prediction for B^B'^ A*^/^ ) is also indicated. We observe 
that in the cancellation region the Higgs- and Z-penguin magnitudes are comparable while 
the box contribution is negligible. This is suggestive of a cancellation between the second 
and third class of diagrams in Fig. [H To observe this cancellation we individually plot the 
absolute values of the form factors Fs^p and 2m^F^ of Eqs. fl2.6H2.9[ ) in the right panel of 
Fig. [3l At the minimum point of the total branching ratio (thick-dashed line in left panel 
of Fig. [3]) \Fp\ is approximately equal to |2m^F^| and \Fs\ is negligibly small. This can be 
explained from the form of Eqs. (12. 6p and (12. 7p . If one assumes S'^^r = (^di/?)*) then Cslr 
and CsRLi the two Wilson coefficients most sensitive to the variation of have similar 
sizes and opposite sign and thus interfere destructively in the amplitude. 

Bounds on the 5 parameters governing squark flavour mixing have been presented in 
the literature using the mass insertion approximation (MIA). In particular, Refs. [38] and 
[39] bound \5'f'j^i\ < 0.3 and ^ 0.02 for a particular point in the parameter space, 

TTiq = M3 = 350 GeV. On the other hand, the results in Fig. [3] arise from an extensive 
scan of the experimentally allowed parameter space without resorting to MIAo. Thus the 



•^Note that references to the i5-parameter in this paper are mainly for comparison and presentation. 
Any other parameter that characterizes the squark mixing would also be appropriate. Recall that our 
calculation is not based on expanding this parameter around zero and keeping only leading terms (MIA 
approximation). Instead, we numerically diagonalize all relevant squark matrices and plug the result into 
the expressions given in the Appendix. 
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bounds on the 5s presented here are both different and more representative of the range 
of possibihties in the general MSSM. The results of this scan show that S'^^^^ is still rather 
weakly constrained, whereas ^Jl^ < 0.08. 

We remark here that varying or has almost no effect on B{Bg fJ'^f^^) which 
takes values along a narrow band. 

3.4 Predictions for B{Bd A^^A^ ) 

We present the corresponding MSSM predictions for B{Bci Ai"*"/^") in Fig.lHwhere 5^"^^ or 
^dL_R varied instead of S'^^j^ or S'^\j^ along with the other SUSY parameters in Table O 
Some sequences of points disappear due to the experimental constraints given in Table [3l 
Note that varying ^^l^.; ^dLR almost no effect on B{B^ — > ji'^jj,^). 

For both cases there exist points where B{Bd f^^f^~) is reduced by an order of 
magnitude relative to the SM. These points are more sensitive to low tan/5 in the 'LL' 
case and fall into the case of Eq. ( 13. 3p . It is also interesting to look at the the ratio 
B{Bd fi^ fi^)/B{Bs yU^/^ ) versus and 6li\j^ , plotted in the right panel of Fig. HI 
Unlike the Standard Model where B{Bd fi+^-)/B{Bs n^fi') ^ \Vtd/Vts\^ < 0.03, 
the MSSM can enhance this ratio by a factor of ten even for small values of 5^^^ or 
^Flr- This suggests that collider searches for B{Bd /i^yU^) are as important as those 
for B{Bs — > fi^fi^). This observation has been already discussed in the literature [24] 
in the leading tan (3 approximation. On the other hand MSSM can further reduce the 
ratio in the 'LL' case by an order of magnitude due to the aforementioned cancellations in 

4 Conclusions 

We have presented a complete, one-loop calculation of the branching ratios for the rare 
decay modes — without resorting to the limits of large tan/3, minimal flavour 
violation or SUSY breaking scale dominance. Our final expressions are presented in an 
appendix and are also available as a computer code (see footnote 2). We have used this 
code to perform a numerical exploration of the MSSM parameter space for the modes 
B^^ —>■ fi'^fi~. We find that there exist cancellation regions where the contribution of 
diagrams with supersymmetric intermediate particles interferes destructively with purely 
Standard Model diagrams, thus allowing the branching ratio to be significantly smaller than 
the Standard Model prediction. We identify possible mechanisms of such cancellations and 
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Figure 4: Upper panel: Predictions for B{Bd versus tan/5 (lei^) and B{Bd —> 

fi~^fi~)/B{Bs f^^f^~) versus 5^j^^(right) from the scan of MSSM parameters in Table [2] 
and under the constraints of Table [31 The solid line shows the upper bound from the 
Tevatron and the dashed line shows the SM expectation. Lower panel: Similar to the 
upper panel but with varied. In all plots only the 6 indicated is varied with all other 
mass insertions set to zero. 



explain why they can occur for certain regions of parameter space. If supersymmetry is a 
proper description of elementary interactions, such effects may effectively hide the dimuon 
decay mode from the LHCb even though it is supposed to be one of the experiment's 
benchmark modes. We have also shown that, barring the cancellations mentioned above, 
supersymmetric contributions in the general MSSM typically tend to enhance the branching 
ratio for B^ ^ —* ^^fi~ even for moderate values of tan/? < 10 so that an experimental 
measurement close to the SM prediction would put strong bounds on the size of allowed 
flavour violation in the squark sector. Finally, we show that the B^ fi^fi^ decay can also 
be either suppressed or enhanced compared to its SM expectation, leading in some cases 
to a situation where the rate of the B^ decay is larger then that of the B^. 

Our analysis provides a quantitative assessment of the viability of non-minimal su- 
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persymmetric flavour structure and its consequences in the neutral S-meson dimuon decay 
modes in light of existing experimental constraints. This is especially relevant due to recent 
experimental hints for non-minimal flavour structure between the second and third quark 
generations [40,41]. Further, there have also been recent model-building analyses of super- 
symmetric models not constrained to the "Minimal Flavour Violation" scenario [42-44]. 

We conclude that new physics, in particular non-minimal flavour violating supersym- 
metry, can manifest itself at future experiments as either an enhancement or a suppression 
of the Bg — > decay rate relative to the Standard Model. 
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Appendix A Wilson coefficients 



In this appendix we provide explicit results for the contributions to the ^ ^ £^£~ Wilson coefBcients 
coming from the self-energy, Higgs- and Z-penguin and box diagrams. In general, Wilson coefficients 
defined in Eq. (|2.1[) can be decomposed as: 
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In the expression above Bx arc the box diagram contributions, Fzx are the Z-penguin irreducible (triangle 
diagram) contributions, Fhx and Fax are respectively the irreducible scalar and pseudoscalar Higgs 
penguin contributions and T,dx are the self energy contributions. Indices are assigned as follows: /, J are 
the generation indices of quarks involved in the process, e.g. I,J— (3, 2) for Bg decay and (3, 1) for B^i 
decay, and K, L are the indices of outgoing leptons, e.g. K,L — (2, 2) for B — > fi^fi~ decay etc. For 
the definition of Higgs mixing matrices Zh, Zr, Higgs boson masses and other symbols we refer reader to 



16 



Ref. [45] , the notation of which we use consistently in this appendix. 



Appendix A.l Loop Integrals 

Here we collect the analytic forms of the relevant loop integrals for this work. The two-point loop integral 
Bi is defined as: 



1 



(471)2-^''"''"''"" '^'^ ' J (27r)2 (P-TO2)[(;._^p)2_^2]• 

The explicit formula for the 2-point loop integral Bi at vanishing external momentum is: 

Bi{0,x,y) = 2 + ^C2ix,y,y), 



(A.ll) 



(A.12) 



where C2{x, y, y) is given in eq. (|A.17[) . 

The 3- and 4-point loop integrals at vanishing external momenta are defined as: 
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The explicit formulae are listed below (we give also expressions for some 3-point functions proportional to 
higher momenta powers, useful in Higgs-penguin calculations): 
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where the divergent piece A = + log(47r)7£; — 1 and 11 is the renormalisation scale. 
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Appendix A. 2 Feynman Rules 



We use the following generic Feynman rules for the calculations below (V^, Sa and are generic vector 
bosons, scalars and fermions, respectively): 



Sa 



^{Vfs7P^ + yfs7PR 



' Ra 



Sb 



I 

t 

I 



Rb 



A^rabc 



Explicit formulae for the generic couplings can be inserted from [45]. In our calculations, V can be Z or 
W bozon. The indices Q, i?, S can denote CP-even or CP— odd Higgs bosons (-fff, A^), squarks (C/^, Di) or 
sleptons {Li^v^). The indices /, g can denote quarks {d^,u^), leptons [l^ ,v^) or charginos and neutralinos 



Appendix A. 3 Box Diagram Contribution 

Box contributions to the Wilson coefficients are denoted by P^xy'- ^ labels the operator type, Z ~ S. V.T 
for scalar, vector or tensor, respectively. X and Y label the handedness, X, F e {L,P}. /, J and K,L 
are quark and lepton generation indices, as described at the beginning of the appendix. Here and in the 
following sections, we strictly follow the notation of [45], where expressions for all mixing matrices, vertices 
and other symbols used can be found. 
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Appendix A. 4 Z-penguins 



Fzx are one- loop triangle-diagram contributions to the X-handed {X = L,R) d d coupling. The 
expression below is valid only for the flavour violating case I ^ J since, in order to simplify the formulae, 
we have dropped some terms appearing only for I = J. 



+ 



+ 



+ 



+ 



M=l 



1- 



w 



cw 



J2 K^^K^^^Z]!Yfm^MM'y„M'^,mlj 



M=l 
3 2 



2swcw 



M=l 1=1 



"W 



^ (c2{ml+,ml^,mU - ^ + (l - ^) rnl^C!o{ml+,ml^,mlj 



1=1 m,n=l 



LIlmirLJln* 
dUC ^dUC 



- '^rnc^rnc„V£ScCo{mf;^,mt.^,m'c^ 



/,n=l m=l 

1 ^ ^ 

oE E K 



-frLIlm-f/LJnm* 



ZUU^dUC ^dUC 



LIl'm-\rLJlri-k 
dDN ^dDN 



v.. 



ZNN 



1=1 m,n=l 



/-< / 2 2 2 \ 



1 ^ / 1 

l,n=l m=l ^ 
l,n=l ^ 



(A.32) 



21 



2s 



W 2 



M=l 1=1 
,2 2 



1 - 



4s 



w 



I ( I I 1 \ ( I I ( I I I \ 



rRIhnjrRJln* 



1 

2 ^ ^ ^dUC^'^dUC 
m.n— 1 Z^l 

^ \ ^ \ ^ T/ni rrRIlmrrRJnm* I ri ('™2 _2 _2 \ , \ 

2 ^zuu^duc Vduc [(^2{mc„^,mu^,muJ + - j 

m=l/,n=l ^ ^ 



1 ^ 



RIlra^rRJln-k 
dDN ^dDN 



7n,n— 1 /^l 

1 rt'x/'Lnm* / 2 2 2 \1 



4 6 

2^ 2^ ^ZDD^dDN ^dDN 



m— 1 l.n—1 



E^. 



/ri 7(/+3)/^(J+3)ra* 



ZDD^D 



D 



(^C2{ml,,m%^,ml,J + 



(A.33) 



Appendix A. 5 Higgs penguins 

Fhx and Fax denote the CP-even and CP-odd one-loop triangle-diagram contributions to the X-handed 
{X = L,R) couplings d'cl-^H^ and d^d-^A" {H^ = H°,H^ = h°,Al = A",Al = G°). Appropriate 
expressions are listed below - please note that the explicit factor of "i" in the CP-odd higgs form factors is 
superficial and comes from the definition of vertices in [Appendix A. 2 For the CP-odd Higgs, the relevant 
vertices defined in this way are, for real Lagrangian parameters, purely imaginary so that iV is a real 
number. 
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Appendix A. 6 d self-energy contributions 

Finally we list the formulae for the one-loop down quark self energy contributions: 



2 



3 



2 



^di = ^Y.K"K"^B,{0,ml^,M^ 



L=l k=l 
2 6 

1=1 fc=l 

4 6 

1=1 fe=l 

+ ^E(^^'^i^'*+4'^'^'^l/^'''*)^i(0'™G,mLj (A.38) 



3 .=1 
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2 3 



L=l k=l 

i=\ k=i 

4 6 

^ Sr^Sr^f^rLIkl^rLJkU ^rRIklrrRJkl*\r) (n. 2 2 \ 

;=i fe=i 

/I 2 6 

^9s X^fylkyJk* y{I+3)k ^{J+3)k*.„ ,„ 2 ^2 \ /a qQ\ 



fe=l 



L=l k=l 

2 6 

^ X^'ST' f-t/LIkl-irRJkl* j_ TrRIkl-irLJkl*-^^ r> fn ^2 2 n 

~ 2 ^ '^'^'^ '^'^'^ ~^^duc^duc )TT^CiBQ[0,m(j^.,mjjJ 
1=1 k=i 

1 * ^ 
i=l fc=l 

+ + 4^+^)''^Zi^'=*)B„(0, ml, mj,^ ) (A.40) 



fe=i 



= l{Yj-Y/)j2T.^HZh'm^,Y^^K"K^''*Bo{0,ml^,ml^] 



2 

L=l k=l 

2 6 

2 



1=1 k=i 
1 ^ 

^ V^/-T/L/feiTAflJfci* T.rRIkUrLJkl*\^ r, (n 2 2 \ 

Z^Z^i ^dDN UdN - ^dDN ^dDN )^Ni Bq (U, m^^ , m^,^ ) 



2 

i=l fc=l 

^2 



6 



/c=l 
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